Theory of the Metal-Insulator Transition in PrRu 4 Pi 2 and PrFe 4 Pi 2 
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All symmetry allowed couplings between the 4/ 2 -electron ground state doublet of trivalent 
praseodymium in PrRu4Pi2 and PrFe4Pi2 and displacements of the phosphorus, iron or ruthenium 
ions are considered. Two types of displacements can change the crystal lattice from body-centred cu- 
bic to simple orthorhombic or to simple cubic. The first type lowers the point group symmetry from 
tetrahedral to orthorhombic, while the second type leaves it unchanged, with corresponding space 
group reductions Im3 — >Pmmm and Im3 — »Pm3 respectively. In former case, the lower point-group 
symmetry splits the degeneracy of the 4/ 2 doublet into states with opposite quadrupole moment, 
which then leads to anti-quadrupolar ordering, as in PrFe4Pi2- Either kind of displacement may 
conspire with nesting of the Fermi surface to cause the metal-insulator or partial metal-insulator 
transition observed in PrFe4Pi2and P1-RU4P12. We investigate this scenario using band-structure 
calculations, and it is found that displacements of the phosphorus ions in PrRu4Pi2(with space 
group reduction Im3 -^Pm3) open a gap everywhere on the Fermi surface. 

PACS numbers: 71.27. +a 61.50.Ks 71. 15. Mb 71.70.Ej 



I. INTRODUCTION 

The rare earth-filled skutterudites (RT 4 Pni2; R=rare 
earth; T = Fe, Ru or Os; Pn = P, As or Sb) exhibit 
a wide range of phenomena, including superconductiv- 
ity in nearly all of the La-compoundsi* 2 * 3 -, as well as 
PrRu4Asi2 2 i and PrRu 4 Sbi2 3 >, heavy fermion supercon- 
ductivity in PrOs4Sbi2^, semiconducting or scmimetal- 
lic behaviour in the Ce-compounds 5 and ferromagnetism 
in several of the Nd- and Eu-compoundsi 3 ^*^ and 
UFe 4 Pi2- PrRu 4 Pi2 under goes a metal-insulator (M- 
I) transition at approximately Tmi = 60K, 8 * 9 - while in 
PrFe 4 Pi2 there is a partial M-I transition at a much lower 
temperature, seen as a sharp upturn in the resistivity at 
Tmi = 6.7K jiSiii which then tends to zero at low temper- 
atures. In both cases, the M-I transition is accompanied 
by a structural phase transition which doubles the vol- 
ume of the primitive cellM^ 

In PrFe 4 Pi2, a peak in the magnetic susceptibility 
is observed at Tmi^ and Tmi is field-dependent, but 
there is no magnetic ordering, 14 However, in-field neutron 
diffraction experiments found direct evidence of anti- 
quadrupolar ordering below Tmi^ In contrast, Tmi of 
PrRu 4 Pi2 is not field dependent (as seen in the specific 
heat 1 ^ and thermal expansion-i) and no peak is observed 
in the susceptibility either^ which suggests that there is 
no ordering of any kind at Tmi- An upturn in the sus- 
ceptibility at very low temperatures hints that there may 
be some kind of ordering below T — 0.35K, 1 and is likely 
related to the low-temperature magnetic field dependent 
peaks found in the thermal expansion and specific heat, 



whose positions decrease in temperature with decreasing 
fields 

Nesting of the Fermi surface may be the common fea- 
ture of the M-I transition in both of these materials. 
Band-structure calculations for LaFe 4 Pi2r^ PrFe 4 Pi2 
and PrRu 4 Pi2^ found an approximately cubed-shaped 
hole-like Fermi surface, with nesting wavevector q = 
(1,0,0) (due to the 48th band in LaFe 4 Pi2 and the 
49th band in PrFe 4 Pi 2 and PrRu 4 Pi 2 ). The band itself 
has a roughly flat dispersion, which means that slight 
differences can cause substantially different Fermi sur- 
face topologies. Thus, in PrRu 4 Sbi2 and LaRu 4 Sbi2 
the resulting FS does not have the nesting property 19 : 20 
LaFe 4 Pi2 and PrFe 4 Pi2 possess an additional smaller 
hole-like spherical Fermi surface, so that the M-I tran- 
sition in PrFe 4 Pi2 is incomplete, and it remains a metal 
at low temperatures. It has been suggested that a nested 
Fermi surface is also a possible cause of the resistivity 
upturns seen in NdFe 4 Pi2, SmRu 4 Pi2, GdRu 4 Pi2 and 
TbRu 4 Pi2-2i 

If nesting is a prerequisite for the M-I transition, then 
so is a cell-doubling structural instability, which halves 
the Brillouin zone (BZ) at the nesting wave vector. All of 
the rare earth-filled skutterudites crystallise in the body- 
centred cubic (bec) lattice with space group Im3 (#204), 
except for the low-temperature phases of PrFe 4 Pi2 and 
PrRu 4 Pi2. In these cases, superlattice reflections are a 
clear indication of doubling of the primitive cellii 2 ^ The 
change in structure appears to be due to a displacement 
of Fe ions in PrFe 4 Pi2 and both of Ru and P ions in 
PrRu 4 Pi2, with space groups Pmmm (#47) and Pm3 
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(#200) respectively. 

In Section II, all symmetry allowed couplings of Pr 3+ 
4/ 2 electrons to the lattice are considered. It is shown 
that anti-quadrupolar ordering of the 4/ 2 electrons fol- 
lows when the point group symmetry at the Pr site is 
lowered from tetrahedral to orthorhombic (as is the case 
in the Im3 — > Pmmm transition). In Section III we 
present band structure calculations which show that a 
gap opens everywhere at the Fermi energy when P ion 
displacements (with space group reduction Im3 — > Pm3) 
are considered in PrR^P^i^ We summarise our results 
in Section IV. 



II. STRUCTURAL PHASE TRANSITION AND 
ANTI-QUADRUPOLAR ORDERING 

According to XANES studies on PrR^Pi^ and sus- 
ceptibility measurements on PrFe4Pi2r^ the valency of 
the Pr ion is +3, ie, a 4/ 2 configuration with total an- 
gular momentum J — A. The Ru/Fe ions form a crys- 
tal field with Oh symmetry with respect to the Pr ions, 
but the P ions lower the symmetry to T/,. Under Oh 
symmetry, the J = 4 state splits into a singlet, a dou- 
blet and two triplets, corresponding to the representa- 
tions Ti © T 3 © T 4 © T 5 of O respectively. Under the 
actual symmetry T/j, the r 3 doublet is split into states 
related by time reversal symmetry (complex conjugates). 
Crystal field splitting alone cannot determine the ground 
state, but specific heat measurements on PrFe4Pi2^ and 



PrRu4Pi2 1 ^ favour the doublet, which is 
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These states carry quadrupole moment 



J 2 |r±) =±8h 2 . 
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On the other hand, a singlet ground state, 



|r x ) = v/5724|4) - y/7/U\0) + jhjn - 4) (3) 

with a low lying triplet first excited state has not been 
conclusively ruled out. The exact form of the triplets is 
unknown in T h , because in general they are linear com- 
binations of 
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In the following discussion, we often refer to the repre- 
sentations of Oh, but since the actual crystal symmetry is 
Th, we should use the following correspondences between 
the representations of O and T: 



Dimension 





T 


1 


Ti 


r a 


1 


r 2 


r a 


2 


r 3 


r b ©r c 


3 


r 4 




3 


r 5 


T d 



We use the representations of Oh for convenience, espe- 
cially in the discussion of the displacement modes of the 
Fe/Ru ions, since by themselves they form a crystal field 
with Oh symmetry. However, the above table should al- 
ways be used to reduce the representations of Oh to those 
of T h . 

Now we consider all possible modes of the displace- 
ments. There are eight Fe/Ru ions within the bcc con- 
ventional cell, which are found half-way between neigh- 
bouring Pr ions (the centres of the hexagonal faces of the 
bcc Wigner-Seitz (primative) cell). They form corner- 
sharing cubes about each Pr ion. There are 24 displace- 
ment modes, which transform according to the represen- 
tations riffir^ffir 3 ffir 3 ffir4ffi2r^ffi2r 5 ffir^ ofOh., where 

the primes indicate the odd (ungerade) representations. 

Twelve P ions are found inside the bcc primitive cell 
(note that the bcc primitive cell is half the volume of 
the bcc conventional cell). The 36 modes of the P ion 
displacements transform according to the representations 



5r' d of TV 



Each 
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of the 36 modes on a single site can be either in-phase or 
anti-phase with the modes on the nearest neighbours, for 
a total of 72 modes in the bcc conventional cell. However, 
only the anti-phase modes will cause a cell-doubling; the 
in-phase modes may lower the point group symmetry but 
the lattice structure will remain body-centred. 

Next, we consider all symmetry allowed couplings be- 
tween the 4/ 2 angular momentum and the lattice. The 
angular momentum cannot couple to the odd representa- 
tions, and since pf transforms as r 3 ©T 3 = ri©r 2 ©r 3 
we find that only the Fe/Ru lattice distortions which 
can couple to the 4/ 2 T 3 doublet are the Ti and T 3 
modes. The T± mode is a dilatation of the cube with 
the eight Fe/Ru ions on each corner and is represented 
by the displacement vector v = (1,1,1). v is short- 
hand notation for the 24-component vector of the po- 
sitions for each of the eight ions. The Ti mode does 
not lower the point group symmetry, but it does dou- 
ble the primitive cell and change the lattice from bcc to 
simple cubic (sc). The two-dimensional T 3 mode is rep- 
resented by the vectors v = (1, 1, —2) and v = (1, —1, 0). 
Again, this is shorthand notation for the distortion of 
the cube with eight Fe/Ru ions on the corners into a 
solid rectangle. Such a distortion also changes the lat- 
tice, and lowers the point group symmetry from to 
T>ih or Y>2h respectively. Under crystal field symme- 
try, there is only one coupling for each of the T\ and 
r 3 modes, since the identity (Fi) only appears once in 
the decompositions Ti © T 3 © T 3 = Ti © T 2 © T 3 and 
T 3 © T 3 © T 3 = Ti © r 2 © 3r 3 . Under the actual sym- 
metry T/j, r 2 reduces to r a , and there are two couplings 



3 



for each mode instead of one. 

There are more choices of coupling to the P ion dis- 
placements. There are two different r a modes, repre- 
sented by v — (1, 1, 1), and each of these has two cou- 
plings to the 4/ 2 doublet. Also there are two kinds 
of Tb T c modes, represented by v — (1,1,-2) and 
v = (1,-1,0), which altogether have four different cou- 
plings to the 4/ 2 doublet. 

Although any of the above couplings might lead to a 
M-I transition via coupling to the conduction electrons, 
there is a significantly different consequence for those dis- 
tortions which lower the point group symmetry (namely 
the T 3 or Tb © r c modes) compared to those which do not 
(the Ti or T a modes). When the point group symmetry 
is lowered by the T3 or Tb © T c modes, the degeneracy of 
the 4/ 2 doublet is lifted, and as we shall show next, the 
result is anti-quadrupolar ordering. 

We consider first the 1^ or r a ©r{, modes. The general 
form of the Hamiltonian which describes the coupling 
between Fe/Ru ion or P ion displacements and the 4/ 2 
doublet in the bcc conventional cell is 

h = ec/v+rv* +/' t /'+r t r) 

+\wr a* - f^'*A*) + h.c] + ||ai 2 

+ ||A| 4 + ^L4| 6 + g(A 6 ±^) 

+^(A 6 -A* 6 ). (10) 

The operators p and f'< act on the two neighbouring 
Pr sites and correspond to the complex wavefunctions 

I/) = (|r+) - i|r 3 -»A/2 and 1/*) = (|r+) + i\rj ))/V2. 

A is the complex amplitude of the displacement v — V3 + 
XV2 +X 2 vi of either eight Fe/Ru ions or 24 P ions, where 
X = cxp (^if-) and the displacement vectors are V3 = 
(-1,-1,2), v 2 = (-1,2,-1) and Vi = (2,-1,-1). The 
actual displacement is then 2KeAv. A corresponds to the 
anti-phase mode as described above, therefore it changes 
sign under the bcc lattice translation t = (§ , § , § ), hence 
odd order in A terms do not appear in H. (Note that 
under the same translation / — * /'). All of the coupling 
constants are real, except for T, which is real in Oft, 
but complex in Tft (therefore, there are twice as many 
coupling constants under T/j, as noted above). The C3 
term is prohibited in Oft, but allowed in T/j. C\ should 
be greater than C2 and C3 for stability. 

We begin by diagonalising H with respect to /. Writ- 
ing A — ae la and T — te lT , the eigenvalues are e ± ta 
with eigenvectors e l ^ a+r ^\f) ± |/*). Then the ground 
state projection of H is 

H e rr = —ta+—a +—a + — a + — a cos6a+ — a sin6a 
2 4 6 6 6 

(11) 

We will assume for the moment that r = and C3 = 
(as in O/j symmetry). Then, for C 2 < 0, H e jf has 
minima at a = 0, 7r/3, 2ir/3, it, 47r/3, 5ir/3. The corre- 
sponding displacements are proportional to ±(—1, —1, 2), 



±(— 1, 2, —1) and ±(2, —1, —1). These displacements re- 
duce the crystal field symmetry from to D^, thereby 
splitting the 4/ 2 doublet, so that now the ground state 
is non-degenerate. For example, for the (—1,-1,2) dis- 
placement (a — 0), the 4/ 2 ground state is |r^~) . The 
4/ 2 ground state on the neighbouring site is orthogonal 
to this (because of the change of sign on T), therefore 
the 4/ 2 ground state alternates between |r^") and \T^} 
from site to site. Since the states \Tf ) carry opposite 
quadrupole moments, it follows that the ground state will 
have anti-quadrupolar ordering. The displacement dou- 
bles the primitive cell and changes the structure from bcc 
to simple orthorhombic (Im3 to Pmmm). 

When C'2 > 0, H e ff has minima at a = 
7r/6, 7t/2, 57r/6, 77r/6, 37r/2, ll7r/6. This leads to domains 
with displacements of the form ±(1,-1,0), ±(1,0,-1) 
and ±(0,1,-1) The corresponding 4/ 2 ground states 
carry quadrupole moments J 2 — J 2 etc. which alternate 
sign from site to site in the same way as for C2 < 0. The 
space group reduction is also the same: Im3 to Pmmm. 

When the coupling constant T is complex and C3 is 
non-zero (as they are in general in T^ symmetry), the 
angle r is non-zero and the angle a is no longer fixed. 
The ground state can be a mixture of \T$) and |r^~) and 
the displacements are a mixture of (2, 1, 1) and (1, —1, 0) 
modes. However, for the case when the 4/ 2 doublet is 
more strongly coupled to the Fe/Ru displacements, the 
Oft, symmetry is broken only by the presence of the P ions, 
which interact indirectly through their contribution to 
the crystal field. Therefore, in this case, it may be appro- 
priate to take the Oft — > T^ symmetry lowering as a per- 
turbative effect, and then the resulting anti-quadrupolc 
moment should be close to the maximum value which is 
found under Oft symmetry. 

Now we consider the Ti or r a modes. The Hamiltonian 

is 

H = e(/ t / + /' t /' + /* t /* + /'* t /'*) 
+TA(pf- f'\f + h.c.) +lu A 2 

+ !^ 4 + ^\ (12) 

where A is now the operator for the displacement v = 
(1,1,1) and all of the coupling constants are real. In 
this case, the 4/ 2 ground state remains degenerate and 
there is no lattice distortion. Therefore, in order for this 
type of distortion to occur, the source of the instability 
must be something else, such as coupling to conduction 
electrons, as we show in Section III. 

The assumption that the doublet is the ground state of 
the Pr 3+ ion is not a strict requirement. In the alternate 
scenario of a singlet ground state and low-lying triplet ex- 
cited state, Kiss and Fazekas 2 - showed how to couple the 
r 3 lattice distortion mode to the quadrupole moment of 
the triplet. In this case, the Hamiltonian takes the same 
form as in l|10|l except that the operator / corresponds to 
the wavefunction |/) = (|r|) + xlH) + X 2 \^ »/V3. Di- 
agonalising H breaks the degeneracy of the triplet; the 
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energy of one state is lowered by TA, and another is 
raised by the same amount. If TA is large enough then 
the state e 4( - Q+T) |/) - |/*) is the ground state of the 
system. Note that this state carries no magnetic mo- 
ment, but in general has a quadrupole moment. Anti- 
quadrupolar ordering follows as before. 

Thus, beginning with all possible couplings between 
the Pr 3+ 4/ 2 doublet and the lattice, we find that there 
are only two types of lattice distortions that are consis- 
tent with a cell-doubling structural phase transition. In 
the first case (r 3 ), the point group symmetry is lowered 
from T/j to T)2h, and the doublet is split into states car- 
rying opposite quadrupole moment. Anti-quadrupolar 
ordering is a consequence of the lower point group sym- 
metry and doubling of the primitive cell. For the sec- 
ond kind of cell-doubling distortion (Ti), the point group 
symmetry remains T/j, and there is no quadrupole or- 
dering. Experimental evidence suggests that the first 
scenario applies to PrFe4Pi2 and while the second one 
applies to PrRu4Pi2- 



TABLE I: The calculated band gaps Ae in eV with distor- 
tions. The minus sign means that the valence and conduc- 
tion bands are overlapped. The asterix indicates that a self- 
consistent solution was not obtained. 
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III. BAND STRUCTURE CALCULATIONS 

Band structure calculations were performed using the 
FLAPW-LDA+U method^ to investigate the origin of 
the M-I transition in PrRu 4 Pi 2 . The LDA+U method 
has described the non-magnetic localized / electron 
system^ 

Various Ti-type lattice distortions involving P ion dis- 
placements were considered. The space group is lowered 
from Im3 to Pm3 as the primitive cell is doubled. In 
Im3, we used the observed lattice parameters, i.e. the 
lattice constant a = 8.0424A and the internal param- 
eters u = 0.3576 and v = 0.1444 for P positions. 27 
In Pm3, the positions of Ru ions are equivalent as 8i, 
(1/2 + 8 T , 1/2 + S T , 1/2 + 8 T ). While the P ions occupy 
the inequivalent positions 12j, (0, u + S^, v + 8^) and 12fc, 
(1/2, 1/2 + u + 5%, 1/2 + v + «5f ). When 8 T = 0, 8% = 5* 
and 8^ = 8^, the space group has the higher symmetry 
Im3 with the smaller primitive cell. The site of Pr ions is 
split to two sites (la (0,0,0) and lb (1/2,1/2,1/2)) in Pm3 
due to the inequivalent surroundings. However, since the 
local symmetry m3 for the both sites is unchanged, any 
kind of lifting the degeneracy of the localized 4 f 2 state is 
not expected in Pm3. Here the singlet Pi states for 4 f 2 
electrons as the starting state were assumed and found 
after the self-consistent steps. 

Although Ru ions are experimentally distorted to min- 
imize the total energy, the position of Ru ions does not 
affect the Fermi surface, because the Fermi surface of 
PrRu4Pi2 consists of mainly P-p band. Therefore 8t = 
is assumed in the calculations. Moreover, the anti-phase 
distortion of P ions are considered. Several values of the 
two parameters S u — 8^ = —8^ and 8 V = 8„ = — <5,f 
with 8 U S V < were considered in the calculations, and 
are listed in Table HJ 

Optimal results were obtained for 8 U — 0.003a and 



PrRu 4 P 12 (Pm3) 
«5u=0 . 00a,5v=-0 .00bS) 




X Z M Z X A F A RS X A r £ M T R 

FIG. 1: Insulating band structure for P1-RU4P12 obtained with 
&u = 0.003a and 5 V = -0.004a. 



S v = —0.004a. This distortion produced a gap across the 
entire Fermi surface, as shown in Fig.GJ Fig.|2shows the 
density of states in the vicinity of the Fermi level both 
for Pm3 and Im3 corresponding to Fig^ 
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FIG. 2: The density of states for PrRu4Pi2 for (a) Im3 and 
(b) Pm3 with S u = 0.003a and 8 V = ™0.004a. Solid line, 
dotted line and dashed line indicates total, Ru-d component 
and P-p component, respectively. Note that there are two P-p 
components in (b). 



IV. DISCUSSION AND SUMMARY 

There are significant differences in the nature of the 
transition between PrFe4Pi2 and PrRu4Pi2- The M- 
I transition temperatures differ by an order of magni- 
tude, with T M i = 6.7K for PrFe 4 Pi2±i versus T M i = 
60K for PrRu 4 Pi2 8 . In addition, T M i of the Fe- 
compound decreases in temperature with increasing mag- 
netic fields^SS, while the Ru-compound shows no mag- 
netic field dependence of Tmi (seen as a jump in the 
specific heat)A&, and no anomaly in the magnetic suscep- 
tibility at Tmi® either. This suggests that, in spite of the 
similarities in the shape of the Fermi surface and in the 
structural phase transition, the lattice instability may be 
of a quite different character between the two materials. 

The observation of anti-quadrupolar ordering below 
Tmi in PrFe4Pi2 is strong evidence that the structural 
phase transition in PrFe4Pi2 involves lattice distortion 



modes which lower the point group symmetry from T^ to 
^2h (the modes). Moreover, the experiments indicate 
that it is the Fe ions which are displaced in this mode^ 
As for PrRu4Pi2, the absence of magnetic anomalies 
mentioned at the beginning of this section, as well as 
Raman studies^ and electron diffractioniS, indicate that 
the structural phase transition at Tmi does not change 
the local Pr site symmetry and is not accompanied by 
quadrupole ordering, hence it involves those modes which 
do not lower the point group symmetry (the Ti modes). 
Our band structure calculations favour P ion displace- 
ments over Ru ion displacements since the former are 
shown to be a source of the M-I transition. 

The low temperature upturn in the susceptibility and 
the magnetic field dependence of the thermal expansion 
and specific heat in PrRu 4 Pi2 at low temperatures sug- 
gest that there may be an ordering at a lower tempera- 
ture. It could be due to a second structural phase transi- 
tion which lowers the point group and lifts the degener- 
acy of the doublet, but so far there are no experimental 
clues about the type of deformation involved. There are 
two possibilities, either in-phase or anti-phase modes of 
the T3 type distortions described in detail above. In- 
phase modes were excluded as a source of the M-I tran- 
sition because they do not double the primitive cell. 
Since the cell is already doubled below the first tran- 
sition, exclusion of the in-phase modes is no longer nec- 
essary. Then the space group is lowered through the se- 
quence Im3 — >Pm3 -^Pmmm, but now the 4/ 2 state has 
quadrupolar ordering below the second transition. Only 
the P-ions participate in the in-phase modes, but the 
anti-phase modes may involve either P ions or Fe/Ru 
ions. The anti-phase modes lead to anti-quadrupolar or- 
dering and lower the space group through the same se- 
quence Im3 -^Pm3 -^Pmmm. Given the similarities be- 
tween the transition temperatures and the magnetic field 
dependences of structural phase transition in PrFe4Pi2 
and the proposed second transition in PrRu4Pi2, it seems 
more likely that the second transition in PrRu4Pi2 would 
be due to displacements of the Ru ions, rather than cither 
in-phase or anti-phase modes of P ions. 

To summarise, we have studied structural phase tran- 
sitions coupled to localised 4/ 2 states and to band elec- 
trons in PrFe4Pi2 and PrRii4Pi2. In both cases the 
source of the structural instability is identified, but the 
mechanisms are quite different. In PrFe4Pi2, a Jahn- 
Teller type mechanism occurs, as coupling between P3- 
type displacements and the 4/ 2 electrons lifts the degen- 
eracy and lowers the energy of the ground state of the 
4/ 2 electrons. The space group is reduced from Im3 to 
Pmmm. Coupling to the band electrons removes part 
of the Fermi surface. In PrRu4Pi2, our bandstructure 
calculations show that coupling between lattice displace- 
ments and band electrons will remove the Fermi sur- 
face everywhere and produce a metal-insulator transi- 
tion. The space group is lowered from Im3 to Pm3. This 
produces no change in symmetry at the Pr sites, hence 
no orbital ordering occurs. 
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